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ABSTRACT 


During  the  August  1992  sea  trial  off  the  west  coast  of  Vancouver  Island,  onboard  the 
CFAV  Endeavour,  data  was  collected  from  a  relatively  new  type  of  GPS  (Global  Positioning 
System)  receiver  which  measures  attitude  (heading,  pitch  and  roll)  as  well  as  the  usual  GPS 
measurements  of  position  and  velocity.  GPS  attitude  has  the  potential  to  provide  rapid  and 
precise  attitude  for  a  relatively  low  cost.  Attitude  is  conventionally  determined  using  inertial 
equipment,  which  for  precise  work  is  very  expensive,  large  and  requires  a  long  "settling  time" 
after  tum-on.  The  receiver  being  evaluated  on  this  sea  trial  was  an  Ashtech  Three-Dimensional 
Direction  Finding  System  (3DF).  This  report  provides  some  general  background  on  GPS 
attitude  technology,  presents  the  attitude  errors  observed  during  this  sea  trial  and  characterizes 
these  errors  as  stochastic  processes, 


r£sum£ 


Au  mois  d'aout  1992,  lors  d'un  essai  en  mer  au  large  de  la  cote  ouest  de  Pile  de 
Vancouver,  a  bord  du  CFAV  Endeavour,  des  donnees  ont  ete  mesurees  a  l'aide  d'un  systeme 
relativement  nouveau  de  recepteur  SPG  (Systeme  de  position  globale)  qui  mesure  Pattitude 
(direction,  inclinaison  et  roulement)  en  plus  des  mesures  SPG  normales  de  position  et  de 
velocite.  L'attitude  SPG  peut  permettre  de  donner  rapidement  et  precisement  Pattitude  a  un 
cout  relativement  bas.  L'attitude  est  determinee  conventionnellement  en  utilisant  Pequipement 
inertiel,  qui,  lorsqu’employe  pour  un  travail  de  precision,  s'avere  tres  dispendieux,  est 
encombrant  et  necessite  une  longue  periode  de  rechauffement  apres  Pallumage.  Le  recepteur 
evalue  lors  de  cet  essai  en  mer  etait  un  Systeme  a  trois-dimensions  pour  determiner  la  direction 
de  Ashtech.  Ce  rapport  presente  de  l'information  generate  sur  la  technologie  de  Pattitude  SPG 
ainsi  que  les  erreurs  d'attitude  observees  durant  l'essai  en  mer.  II  caracterise  ces  erreurs  come 
etant  des  procedures  stochastiques. 
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EXECUTIVE  SUMMARY 


During  the  August  1992  sea  trial  off  the  west  coast  of  Vancouver  Island,  onboard  the 
CFAV  Endeavour,  data  was  collected  from  a  new  type  of  GPS  (Global  Positioning  System) 
receiver  which  measures  attitude  (heading,  pitch  and  roll)  as  well  as  the  usual  GPS 
measurements  of  position  and  velocity.  This  type  of  receiver  uses  GPS  carrier  phase  difference 
measurements  between  pairs  of  antennas  (using  interferometric  techniques).  This  method  has 
been  discussed  in  the  literature  since  the  late  1970's  but  is  only  now  becoming  practical  due  to 
recent  developments  in  receiver  hardware  technology 

Attitude  is  conventionally  determined  using  inertial  equipment,  which  for  precise  work 
is  very  expensive,  very  large  and  heavy,  requires  a  long  "settling  time"  after  turn-on  and  has 
high  latitude  limitations.  This  new  type  of  GPS  receiver  however  has  the  potential  to  provide 
rapid  and  precise  attitude  for  a  relatively  low  cost.  Of  course  inertial  attitude  also  has  several 
advantages  over  GPS,  such  as  providing  high  bandwidth  output  (high  data  rates)  and  being 
unaffected  by  jamming,  spoofing,  multipath,  or  loss-of-signal.  It  would  seem  that  the  best 
solution  in  many  cases  would  be  an  integrated  system  with  GPS  attitude  combined  with  a  low 
cost  inertial  system.  To  investigate  the  potential  performance  of  such  an  integrated  system  it  is 
necessary  to  characterize  the  errors  of  the  GPS  attitude  in  a  stochastic  model 

The  receiver  being  evaluated  on  this  sea  trial  was  an  Ashtech  Three-Dimensional 
Direction  Finding  System  (3DF),  which  is  a  24  channel,  C/A  code  receiver  with  four  antennas 
For  this  trial  the  antennas  were  configured  as  a  roughly  level  square  with  10  metre  sides. 
Precise  position  and  attitude  reference  data  was  available  throughout  the  trial  from  the  DIINS 
(Dual  Inertial  Integrated  Navigation  System)  equipment,  which  was  also  onboard  The 
accuracy  of  the  raw  inertial  reference  attitude  was  at  least  about  0.07  degrees  in  heading  and 
0.03  degrees  in  pitch  and  roll  (rms).  More  importantly,  the  noise  and  other  short-correlation¬ 
time  errors  of  the  reference  system  are  much  smaller  than  that  of  the  3DF.  The  reference  data 
was  available  at  0.16  second  intervals  throughout  the  7  day  trial,  providing  almost  four  million 
data  sets.  Post  processing  the  inertial  reference  data  with  differential  GPS  position 
measurements  improved  the  reference  attitude  accuracy  to  about  0.01  degrees  (rms). 

The  GPS  attitude  was  generally  available  at  one  second  intervals,  but  about  1 0%  of  the 
data  was  not  available  and  a  further  10%  was  flagged  as  invalid  (carrier  phase  cycle  ambiguity 
not  resolved)  and  in  fact  had  very  large  errors.  In  this  case  the  3DF  sets  the  pitch  and  roll  to 
zero  and  the  heading  is  based  on  code  phase  measurements  (which  are  much  noisier  than 
carrier  phase  measurements). 

The  accuracy  of  the  valid  GPS  attitude  data  was  first  measured  at  the  discrete  one 
second  pc'nts  for  which  it  was  available  and  was  found  to  be  about  0.5  degrees  in  heading  and 
1 .2  degrees  in  pitch  and  roll  (rms).  However  this  was  due  to  a  small  amount  of  poor  data  (less 
than  one  percent)  and  perhaps  a  better  measure  of  performance  is  the  95  percentile,  which  was 
0.1  degrees  for  heading  and  1.2  degrees  for  pitch  and  roll.  Furthermore,  most  of  these  pitch 
and  roll  errors  were  dominated  by  biases,  which  presumably  could  have  been  removed  by 
proper  calibration.  When  these  biases  were  removed  the  remaining  pitch  and  roil  errors  were 
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significantly  reduced  to  about  0. 1  degrees  (95%).  As  mentioned,  these  numbers  are  based  on 
the  comparison  of  only  the  available  (at  one  second  intervals)  and  valid  GPS  data  with  the 
corresponding  reference  data  and  are  therefore  referred  to  as  "discrete-dynamic"  errors  Thus 
the  discrete-dynamic  attitude  accuracy  in  this  configuration  appears  to  be  about  2  milliradians 
at  the  95%  level. 

The  true  "continuous-dynamic"  performance  is  somewhat  worse  when  consideration  is 
taken  of  the  changes  in  the  platform  attitude  between  the  one  second  GPS  data  points.  This 
performance  is  measured  by  comparing  the  (unextrapolated)  GPS  data  to  all  of  the  0  1 6  second 
reference  data.  The  rms  errors,  again  dominated  by  a  small  amount  of  very  poor  data,  were  at 
the  0.6  degree  level  for  heading  and  1.3  degree  level  for  pitch  and  roll.  The  95  percentile 
continuous-dynamic  accuracy  of  the  GPS  attitude  was  found  to  be  about  0  4  degrees  in 
heading  and  1.9  degrees  in  pitch  and  roll.  With  the  biases  removed,  these  continuous-dynamic 
pitch  and  roll  errors  were  reduced  to  0.9  degrees  (95%).  The  GPS  receiver  under  evaluation  is 
supposed  to  be  capable  of  a  slightly  higher  output  rate  (2  Hz).  Although  higher  rates  may  lead 
to  increased  noise,  they  could  also  presumably  reduce  the  continuous-dynamic  errors,  and  in 
any  case,  if  this  were  used  to  aid  an  inertial  system  the  continuous-dynamic  errors  would  not 
be  relevant. 

The  stochastic  behaviour  of  the  attitude  errors  collected  over  this  week-long  trial  could 
be  modeled  quite  accurately  as  the  sum  of  three  independent  components:  a  constant  bias, 
uncorrelated  white  noise  and  a  first  order  Markov  process.  The  bias  was  on  the  order  of  one 
degree  and  is  likely  an  installation/calibration  error  which  could  be  quite  easily  removed  The 
noise  and  Markov  processes  were  each  at  about  the  0.03  degree  level  (0.5  milliradian),  one 
sigma.  For  applications  where  this  is  not  acceptable,  or  where  better  bandwidth  is  required 
(such  as  for  fire  control)  then  integration  with  even  a  low  cost  inertial  system  could  provide 
very  good  sub-milliradian  attitude  data. 

Considerable  care  must  be  taken  when  extending  these  results  to  other  situations.  In 
the  static  situation,  such  as  artillery  pointing,  the  noise  may  be  reduced  by  longer  sampling 
periods,  however  the  multipath  errors  may  create  worse  temporally  correlated  errors  (i.e. 
Markov).  In  the  marine  environment  these  multipath  errors  are  partially  randomized  or 
"whitened"  by  the  wave-induced  motion  of  the  ship. 

This  type  of  receiver  therefore  shows  considerable  promise  as  an  aid  to  navigation, 
either  by  itself  if  high  bandwidth  and  continuous  data  is  not  essential,  or  integrated  with  an 
inertial  system. 
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INTRODUCTION 


A  sea  trial  was  conducted  in  August  1992  by  the  Navigation  Group  of  the  Electronics 
Division  of  the  Defence  Research  Establishment  Ottawa.  One  of  the  secondary  objectives  of 
this  trial  was  to  evaluate  the  performance  of  a  new  type  of  GPS  (Global  Positioning  System) 
receiver,  which  is  designed  to  measure  attitude  (heading,  pitch  and  roll)  as  well  as  the  usual 
GPS  measurements  of  position  and  velocity.  This  type  of  receiver  can  compute  attitude  by 
taking  double-difference  carrier  phase  measurements  from  several  (at  least  three)  GPS  satellite 
signals  at  several  (at  least  three)  different  antennas.  Since  the  antennas  are  usually  more  than 
one  wavelength  apart,  there  is  a  significant  cycle-ambiguity  problem  to  solve  as  well 
Sometimes  a  fourth  antenna  is  used  to  assist  in  this  ambiguity  resolution.  This  technique  is 
adequately  described  in  the  literature  (see  for  example  references  [1]  to  [10])  and  is 
summarized  in  Section  2  below. 

This  new  GPS  technology  has  many  potential  applications.  It  could  partially  replace, 
or  at  least  complement,  the  very  large  and  costly  marine  inertial  systems  currently  used  to 
measure  attitude  on  ships.  It  could  provide  rapid  and  precise  alignment  of  artillery  pieces.  It 
could  provide  heading  information  at  extremely  high  latitudes  where  inertial  systems  cannot  It 
could  also  be  used  to  precisely  align  sensors  on  airborne  or  spacebome  platforms. 

Platform  attitude  is  required  for  fire  control  purposes  and/or  sensor  alignment  on  most 
military  platforms,  and  in  particular  for  naval  vessels.  This  is  generally  considered  a  mission 
critical  function,  so  that  highly  accurate,  reliable  and  timely  (i.e.  low  latency)  attitude  data  is 
required.  The  inertial  sensors  currently  used  to  provide  this  data  are  very  costly,  massive  and 
power  consuming.  A  marine  INS  (Inertial  Navigation  System)  typically  costs  several  hundred 
thousand  dollars,  weighs  several  hundred  kilograms  and  consumes  hundreds  of  watts  of  power 
For  high  latitude  operations  it  is  important  to  recognize  that  the  heading  accuracy  of  an  inertial 
system  will  be  seriously  degraded,  increasing  roughly  as  secant( latitude),  thus  becoming 
essentially  useless  in  the  immediate  vicinity  of  the  pole.  The  attitude  determining  GPS 
receiver,  by  comparison,  costs  less  than  $80K  (potentially  much  less),  weighs  less  than  1 0  kg, 
consumes  only  about  12  watts  and  is  not  degraded  at  high  latitudes. 

Unfortunately  there  are  also  expected  to  be  shortcomings  with  the  GPS  technology  in 
the  area  of  performance  (accuracy,  availability,  reliability,  data  htency).  Of  course  GPS  relies 
on  external  line-of-sight  signals  and  is  therefore  subject  to  jamming,  spoofing,  shading  and 
multipath.  GPS  attitude  determination  also  requires  the  installation  of  several  antennas  (at 
least  three),  with  associated  cables. 

It  may  be  possible,  however,  to  overcome  most  if  not  all  of  these  shortcomings  by 
properly  integrating  the  GPS  attitude  receiver  with  a  lower  cost  (and  therefore  lower 
performance)  inertial  system  or  AHRS  (attitude  and  heading  reference  system).  The  degree  to 
which  this  will  meet  the  accuracy  performance  requirements  will  depend  upon  the  error 
characteristics  of  the  INS  and  of  the  GPS  attitude,  the  latter  of  which  is  the  subject  of  this 
report.  The  reliability  and  environmental  requirements  would  also  require  special 
consideration,  but  these  are  not  addressed  here. 
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Table  1  below  compares  the  general  characteristics  of  the  two  competing  technologies, 
for  the  purpose  of  attitude  determination  in  the  marine  environment. 

TABLE  1.  GENERAL  CHARACTERISTICS  OF  INS  vs  GPS 


FACTOR 

MARINE  INERTIAL 

GPS 

Cost  ($K) 

>250 

<  80 

Environmental 

full  mil  spec. 

commercial  rugged 

Heading  Accuracy  (degrees  RMS) 

<  0.05  x  secant(lat-) 

mmm 

Pitch,  Roll  Accuracy  (degrees 
RMS) 

<0.025 

TBD  (=0.1) 

Data  Rate 

400  Hz 

1  Hz 

Time  to  First  Fix 

hours 

seconds 

Size  (m^  ) 

>0.3  (0.5  x  0.5  x  1.5) 

=  0.004  (0.2  x  0.2x0  1) 

Weight  (kg) 

>300 

<5 

Power  (watts) 

>500 

12 

Other  Considerations 

self-contained 

cabling,  antennas  with 
clear  view  of  sky, 
multipath,  jamming,  denial 

Comparison  with  ground  based  or  airborne  inertial  systems  will  be  quite  different,  with 
the  INS  being  considerably  smaller,  lighter,  less  expensive,  less  power  consuming,  and 
requiring  less  alignment  time,  but  also  having  lower  accuracy.  The  GPS  error  characteristics 
will  also  be  different  in  the  land  and  air  environment,  since  these  errors  will  depend  largely  on 
the  platform  dynamics,  the  antenna  baseline  lengths  and  the  multipath  environment.  In 
particular  static  accuracy  can  be  significantly  better  than  dynamic  accuracy  (provided  that 
multipath  is  carefully  avoided)  and  longer  baselines  will  produce  better  accuracy  than  shorter 
ones. 

The  purpose  of  this  report  is  therefore  to  produce  an  initial  characterization  of  the 
dynamic  attitude  performance  of  an  Ashtech  Three-Dimensional  Direction  Finding  System 
(3DF)  in  the  marine  environment,  as  well  as  to  provide  some  general  background  on  this 
general  class  of  GPS  receiver  equipment.  The  Ashtech  3DF  is  a  24-channel  single  frequency 
GPS  receiver  with  four  micro  strip  antennas.  References  [6]  and  [10]  provide  a  more  complete 
description. 
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2. 


BACKGROUND 


The  determination  of  heading  and  attitude  from  GPS  has  been  discussed  in  the 
literature  for  some  time.  Written  over  fifteen  years  ago,  reference  [1]  describes  the  general 
principles  behind  attitude  determination  using  phase  comparison  between  pairs  of  antennas 
The  basic  concept  is  as  shown  in  Figure  I . 


FIGURE  1.  INTERFEROMETER  PRINCIPLE 


The  fundamental  interferometer  equation,  derived  from  Figure  1,  relates  the  measured 
phase  difference  <p  at  the  two  antennas  to  the  angle  of  incidence  9  of  the  signal  to  the 
antenna  baseline: 


d  sin9  =  rik  + 


(1) 


where  d  is  the  baseline  length,  X  is  the  signal  wavelength  (where  in  the  case  of  C/A  code 
GPS  X  =  19  cm)  and  n  is  the  integer  number  of  wavelengths  (phase  ambiguity)  needed  to 
ensure  that  measured  phase  difference  tp  is  in  the  range  -it  <  <p  <  7t  .  This  phase  ambiguity 
can  arise  when  d  >  XI2,  which  is  usually  the  case  as  we  shall  see. 


Assuming  that  this  phase  ambiguity  has  been  resolved,  then  from  equation  1  it  can  be 
seen  that  a  phase  measurement  error  of  5<p  will  result  in  an  angular  measurement  error  of  50 
where: 


59 


2ra/cos0 


6(f> 


(2) 


This  shows  how  the  phase  error  6<p,  due  to  either  receiver  phase  measurement  error  or 
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multipath  error,  will  cause  an  error  in  the  angular  measurement.  Thus  the  important  factors  to 
control  when  determining  the  angular  measurement  error  66  are  the  phase  error  6a,  the 
baseline  length  d  (since  X  is  a  constant  and  6  is  determined  by  the  satellite  geometry  which 
cannot  be  controlled)  and  of  course  the  phase  ambiguity  resolution.  Incidentally,  the  satellite 
geometry,  represented  by  the  cos0  term,  is  generally  such  that  horizontal  angular 
measurements  (heading)  are  more  accurate  than  vertical  angular  measurements  (pitch  and  roll) 
by  a  factor  of  1.5  to  2,  just  as  the  horizontal  GPS  position  measurements  are  generally  more 
accurate  than  the  vertical  GPS  position  measurements,  by  the  same  amount 

In  1977  it  was  suggested  in  reference  [1]  that  this  type  of  receiver  could  attain 
accuracies  of  2  to  20  minutes  of  arc  (0.03  to  0.33  degrees)  with  measurement  rates  of  10  per 
second.  This  was  an  analytical  discussion,  in  the  context  of  an  airborne  platform,  assuming  a 
three  antenna  interferometer  (the  minimal  configuration).  The  general  discussion  in  this 
reference  illustrates  how  the  expected  attitude  errors  due  to  phase  measurement  error  and 
multipath  error  (the  two  dominant  sources  of  error)  are  inversely  proportional  to  the  baseline 
length.  This  discussion  assumed  a  GPS  carrier  phase  measurement  accuracy  of  0  6  degrees 
Since  the  wavelength  of  the  LI  carrier  signal  is  19  0  cm,  this  is  equivalent  to  19x0.6/360  = 
0  03  cm.  This  also  assumed  four  metre  baselines  between  three  antennas  The  suggested 
application  of  this  technique  was  to  eliminate  the  long  alignment  times  required  by  inertial 
systems  on  moving  platforms  (i.e.  on  an  aircraft  carrier  or  in  flight).  This  paper  also  foresaw 
the  integration  of  this  type  of  receiver  with  an  inertial  system. 

Reference  [2],  apparently  from  1981,  provides  a  more  complete  discussion  of  error 
sources  and  provides  simulation  resu'*  suggesting  that  an  accuracy  of  50  prad  (0  003  degrees) 
may  be  achievable.  However  this  assumed  a  static  three  antenna  baseline,  no  multipath  error , 
four  GPS  satellites  visible  and  15  minutes  of  Kalman  filtering  with  an  INS  ( 1  nautical  mile  per 
hour  type).  Ignoring  multipath,  especially  in  a  static  situation,  will  result  in  significantly 
underestimated  errors.  The  Kalman  filter  based  integration  with  an  INS  as  briefly  described 
here,  is  an  excellent  idea,  especially  in  the  dynamic  situation.  (The  assumed  baseline  lengths 
were  1  metre  and  the  GPS  carrier  phase  measurement  accuracy  was  assumed  to  be  0  36 
degrees.)  This  paper  also  concludes  that  0.06  degrees  heading  accuracy  can  be  achieved  with 
5  minutes  of  filtering. 

In  1983  DREO  sponsored  a  proof  of  concept  study  to  investigate  the  feasibility  of 
measuring  azimuth  (heading)  using  this  technique  with  the  older  Transit  satellite  system  The 
technique  used,  called  AZTRANS  (AZimuth  from  TRANSit),  was  originally  developed  by  the 
Applied  Physics  Lab  in  1971,  and  was  adapted  by  JMR  Instruments  Ltd  in  1983  Using  a 
JMR-A1  receiver  with  two  antennas  on  a  static  four  metre  baseline,  it  was  shown  in  reference 
[3]  that  one  degree  azimuth  accuracies  (rms)  could  easily  be  obtained,  and  0  5  degree  accuracy 
(rms)  was  feasible  Note  that  since  Transit  operates  at  400  MHz,  its  carrier  wavelength  X  is 
about  four  times  longer  than  that  of  the  GPS  LI  carrier  (which  is  at  1540  MHz),  leading  to 
proportionally  larger  Transit  attitude  errors  58  for  a  given  phase  measurement  error  5<p, 
according  to  equation  2  above 

Later  in  1983,  reference  [4]  provided  another  discussion  of  this  technique  as  applied  to 
GPS,  with  simulations  suggesting  that  2  4  mrad  (0.14  degrees)  rms  attitude  accuiacy  can  be 
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achieved  with  2  metre  baselines,  a  5.0  mm  phase  measurement  accuracy,  and  three  antennas 
It  also  suggests  that  0.6  mrad  (0.03  degrees)  rms  attitude  accuracy  can  be  achieved  with  20 
metre  baselines. 

As  described  in  reference  [5],  an  attitude  GPS  receiver  was  actually  developed  to  the 
ADM  stage  in  1986  by  Stamav  Corporation  of  Ottawa  (with  support  from  the  National 
Research  Council  of  Canada),  which  demonstrated  0.6  degree  attitude  error  standard  deviation 
(plus  some  bias)  in  static  mode  This  used  a  3  antenna  array  with  1  metre  baselines.  The  failure 
of  this  receiver  to  enter  production  was  blamed  on  the  Challenger  disaster  of  1986,  which  set 
back  the  GPS  satellite  launch  schedule,  delaying  market  demand  for  receiver  equipment. 

The  receiver  evaluated  on  this  trial  is  the  Ashtech  3DF.  The  manufacturer  has 
published,  in  reference  [6],  a  paper  describing  the  3DF  system  and  presenting  some 
performance  test  results.  Most  of  the  results  are  from  static  tests  and  indicate  accuracies 
consistent  with  Figure  2.  For  example  the  claimed  static  accuracy  was  as  good  as  0  006 
degrees  (0.1  mrad).  However  this  used  a  baseline  length  of  110  metres  and  apparently  also 
used  heavy  temporal  filtering.  It  therefore  cannot  be  applied  to  the  dynamic  situation  Shorter 
baselines  of  one  metre  produced  accuracies  of  about  0.3  degrees,  again  in  static  mode,  with 
heavy  temporal  filtering.  The  only  dynamic  data  presented  by  the  manufacturer  was  not 
conclusive  because  of  poor  reference  data,  however  it  suggests  an  accuracy  of  between  0  4 
degrees  and  0.05  degrees.  The  text  of  [6]  suggests  the  latter  figure  but  the  data  presented 
does  not  appear  to  support  this  conclusion  and  in  fact  suggests  something  closer  to  0  2 
degrees.  Again  there  seems  to  be  some  temporal  filtering  involved,  which  hides  some  of  the 
true  dynamic  error.  In  fact,  as  mentioned  in  [6],  three  small  Kalman  filters  (2  states  each)  are 
used  to  smooth  the  heading,  pitch  and  roll  measurements  respectively. 

Reference  [7],  from  1991,  describes  an  attitude  determining  system  built  by  Adroit 
Systems  Inc.  in  1990,  designed  for  static  pointing  of  artillery.  This  system  uses  a  modified 
Texas  Instruments  TI-420  receiver  with  three  co-linear  antennas  to  measure  azimuth  and 
elevation  (heading  and  pitch)  only.  The  antennas  form  one  short  (about  14  centimeters)  sub¬ 
wavelength  baseline  for  phase  ambiguity  resolution,  and  one  longer  (one  metre)  baseline  for 
attitude  measuring.  Since  this  was  a  static  implementation,  the  measured  performance  could 
be  improved  by  extending  the  measurement  interval.  The  azimuth  accuracy  varied  from  about 
0.6  degrees  for  a  one  second  sample  to  0.3  degrees  (standard  deviation)  for  30  second 
processing  to  (apparently)  0. 1  degree  for  360  second  measurements.  (Reference  [7]  in  most 
cases  presents  the  average  errors  rather  than  the  rms  errors,  which  is  not  very  useful.)  This 
reference  also  presents  some  results  from  two  competing  receivers  of  the  same  class  (static, 
two  dimensional),  one  from  Magnavox  and  one  from  Texas  Instruments.  Because  of  different 
implementation  decisions  the  error  characteristics  of  the  three  competing  designs  were 
somewhat  different,  for  example  with  the  Magnavox  errors  having  somewhat  larger  biases  but 
smaller  standard  deviations. 

In  reference  [8]  some  1992  experimental  results  are  presented  comparing  an  Ashtech 
3DF  to  a  Litton  LTN90/100  strapdown  inertial  system.  This  dynamic  test  was  performed  on  a 
pickup  truck  over  straight  and  smooth  highways,  so  that  there  were  actually  very  little 
dynamics.  Care  was  also  taken  to  avoid  multipath,  and  time  of  day  was  chosen  to  maximize 
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GPS  satellite  coverage:  five  satellites  with  Position  Dillution  of  Precision  (PDOP)  between  2 
and  4.  Thus  the  limited  data  presented  in  reference  [8]  (about  40  minutes  of  data)  may  have 
produced  optimistic  results  for  this  1  to  2  metre  antenna  baseline  configuration. 

These  results  from  the  literature  are  summarized  in  Table  2  below. 


TABLE  2.  SUMMARY  OF  LITERATURE  RESULTS 


REF. 

#  OF 
ANTEN 
-NAS 

SIMUL¬ 

ATED 

OR 

REAL 

STATIC 

OR 

DYNA¬ 

MIC 

BASELINE 

LENGTH 

(metres) 

PHASE 

MEASUREMENT 

ACCURACY 

(mm) 

SAMPLE 

PERIOD 

(sec.) 

ATTITUDE 
ACCURACY 
(degrees  rms*) 

HEADING  PITCH/ROLL 

1 

3 

S 

S 

4 

0.3 

0.1 

0.03  -  0.33 

2 

3 

S 

S 

1 

0.2 

300 

0.06 

• 

4 

3 

S 

S 

2-20 

5.0 

720 

0.14-0.03 

0.14-0.03 

5 

3 

R 

S 

1 

0.6 

6a 

4 

R 

S 

1 

0.2 

1 

0.14 

0.25 

6b 

n 

R 

S 

110 

0.2 

1 

0.002 

0.004 

6c 

D 

R 

D 

50 

0.2 

1 

0.05-0.2  1 

- 

n 

3 

R 

S 

1 

30 

0.3 

0.6 

8 

D 

R 

D 

1.5 

0.2 

1 

SB 

0.4 

*  Not  all  references  provided  rms  error  statistics,  and  some  that  do  appear  to  have  confused  rms 
with  standard  deviation,  as  discussed  in  the  two  notes  below. 

NOTE  1.  Concerning  the  heading  accuracy  for  reference  6c  in  Table  2:  The  dynamic  attitude 
errors  of  the  3DF  system  presented  in  reference  [6]  (0.4  degrees)  is  a  mean  rather  than  an 
rms  error.  The  rms  of  course  is  generally  larger  than  the  mean  (according  to  the 
Minkowski  inequality).  Reference  [6]  also  suggests  that  most  of  this  mean  error  (0.34 
degrees)  is  due  to  an  error  in  the  calibration  of  the  reference  gyro  rather  than  to  the  3DF 
system,  and  goes  on  to  discuss  the  resulting  error  in  the  mean  (of  0.05  degrees)  as  though 
this  were  the  dynamic  error.  Unfortunately  this  "error  in  mean"  is  not  a  good  measure  of 
performance  since  it  completely  ignores  the  noise.  Furthermore,  when  the  mean  (i.e.  the 
bias)  is  removed  in  this  way,  basically  all  that  remains  is  noise,  which  from  the  plotted 
data  shown  in  reference  [6]  visually  appears  to  be  at  about  the  0.2  degree  rms  level, 
suggesting  an  rms  value  of  about  0.2  degrees.  (Alternatively  this  could  be  quoted  as  the 
standard  deviation  value,  in  which  case  the  mean  is  removed  by  definition.) 

NOTE  2.  Concerning  the  accuracy  results  for  reference  8  in  Table  2:  The  rms  statistics  given 
for  reference  [8]  appear  to  be  standard  deviation  statistics  (i.e.  with  biases  removed), 
because  in  many  cases  the  "rms”  error  statistic  given  is  smaller  in  magnitude  than  the 
corresponding  mean  error  statistic  (which  is  mathematically  impossible). 
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Table  2  suggests  that  performance  can  be  improved  simply  by  extending  the  antenna 
baseline  lengths  and  the  sample  period.  In  fact  the  basic  relationship  of  equation  2  above  has 
led  various  authors,  such  as  in  reference  [9]  and  the  manufacturer  of  the  equipment  under  test 
(in  reference  [6]),  to  publish  graphs  showing  a  direct  relationship  between  baseline  length  and 
accuracy,  as  shown  in  Figure  2  (which  incidentally  excludes  multipath-type  errors,  which  could 
be  very  significant).  Reference  [10]  contains  a  similar  plot  with  an  additional  line  representing 
the  expected  performance  with  multipath,  indicating  errors  which  are  larger,  by  a  factor  of 
about  2.5,  than  the  non-multipath  errors. 

Unfortunately  there  are  often  problems  with  extending  the  baseline  length  on  many 
platforms.  First  there  is  the  limitation  of  the  platform  size,  and  second  there  is  the  limitation  of 
baseline  stability.  The  true  antenna  baseline  lengths  must  be  known  and  preferably  must  be 
very  stable  under  temperature  fluctuations  and  mechanical  vibrations.  Flexing  on  a  large 
platform  can  therefore  cause  difficulties,  although  it  has  been  suggested  that  this  technique  can 
be  used  to  actually  measure  this  flexing,  as  discussed  in  reference  [10] 


FIGURE  2.  ACCURACY  VS  BASELINE  LENGTH 


Another  method  of  improving  accuracy  is  to  extend  the  sample  period,  allowing  many 
different  measurements  to  be  combined  using  some  averaging  or  filtering  technique  to  reduce 
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the  phase  measurement  error.  This  is  only  possible  when  the  attitude  can  be  assumed  constant, 
as  in  the  static  situation,  such  as  for  artillery  alignment,  or  if  inertial-type  aiding  is  available 

The  effectiveness  of  extending  the  sample  period  (using  temporal  filtering)  is  also 
strongly  determined  by  the  degree  to  which  the  phase  errors  are  temporally  uncorrelated  For 
example  uncorrelated  white  noise  can  be  easily  removed  through  filtering  whereas  a  bias 
(totally  correlated)  cannot.  Although  true  biases  can  be  removed  by  an  installation  calibration, 
errors  which  vary  slowly  in  time  cannot  be  removed  by  either  calibration  or  averaging.  In  the 
static  situation,  multipath  can  produce  these  types  of  errors  which  are  significantly  correlated 
in  time,  limiting  the  effectiveness  of  this  filtering 

Reference  [7]  shows  the  expected  relationship  between  sample  period  and  accuracy,  for 
a  one  meter  antenna  baseline,  as  reproduced  here  in  Figure  3. 


Azimuth 


FIGURE  3.  ACCURACY  VS  SAMPLE  PERIOD 
(1  meter  baseline) 
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One  question  which  comes  to  mind  when  reviewing  this  literature  is  the  long  time  from 
concept  to  realization.  The  likely  explanation  for  this  long  delay  between  the  theoretical 
derivation  of  this  technique  (1977)  and  its  practical  implementation  (1992)  is  that  the  extreme 
miniaturization  of  receiver  hardware,  with  receiver-on-a-chip  and  multichannel  chips,  has  only 
recently  made  it  economical  to  build  into  the  same  receiver  the  large  number  of  receiver 
channels  necessary  to  implement  multi-antenna  attitude  determination.  This  delay  could  thus 
be  partially  blamed  on  the  slow  implementation  of  the  GPS  Space  Segment  (i.e.  the  Shuttle 
disaster),  which  delayed  the  consumer  demand. 
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3.  SEA  TRIAL  METHODOLOGY 


The  general  execution  of  this  sea  trial  is  discussed  more  fully  in  reference  [11].  The 
duration  of  the  trial  was  two  weeks  including  installation  and  removal  of  equipment,  from  August 
14  to  28,  1992.  Seven  days  were  spent  at  sea  collecting  data.  The  3DF  data  set  extends  from 
22:09  on  day  230  to  16:19  on  day  237  (Greenwich  mean  time).  There  were  no  indications  of  any 
unusual  conditions  of  the  GPS  space  or  control  segments  at  that  time. 

The  area  of  operation  was  as  illustrated  in  Figure  4,  which  shows  the  ship’s  track  over  the 
entire  trial,  from  the  differential  GPS  reference  data.  The  sea  trial  was  conducted  on  the  CFAV 
Endeavour,  which  is  a  75  meter,  1600  ton  research  vessel  operated  by  the  Defence  Research 
Establishment  Pacific,  in  Esquimalt  B.C. 


-127  -126  -125  -124  -123 

LAT  degrees 

FIGURE  4.  AREA  OF  OPERATION 
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3.1  Equipment 

Since  the  primary  purpose  of  the  trial  was  to  collect  sensor  data  for  DIINS  (Dual 
Inertial  Integrated  Navigation  System)  development  (see  reference  (12]),  the  navigation  sensor 
complement  included  a  relatively  high  performance  marine  inertial  navigation  system,  referred 
to  here  as  INS1.  This  provides  a  continuous  reference  for  heading,  pitch  and  roll  with  an 
expected  accuracy  significantly  better  than  the  expected  GPS  attitude  accuracy,  essentially  no 
noise  or  transients,  and  a  data  rate  about  six  times  faster.  There  was  also  a  lower  quality 
airborne  INS,  referred  to  here  as  INS2,  which  nevertheless  had  highly  predictable  errors  (RLG 
based,  strapdown,  free  inertial). 

There  was  also  a  differential  GPS  receiver  (DGPS)  providing  position  reference 
accurate  to  about  5  metres.  With  the  DIINS  software,  it  is  possible  to  integrate  this  accurate 
position  data  with  either  of  the  INSs  to  provide  further  enhanced  attitude  data.  This  has  the 
added  advantage  of  providing  covariance  information  (from  the  DIINS  Kalman  filter)  to 
indicate  more  precisely  the  expected  accuracy  of  the  integrated  attitude  solutions. 

The  Ashtech  3DF  4-antenna  array  was  installed  at  the  stem  of  the  ship,  as  illustrated  in 
Figure  5.  These  antennas  formed  an  approximately  square  horizontal  array  with  10  metre 
baselines  in  both  the  longitudinal  direction  (heading  and  pitch  measuring)  and  lateral  direction 
(heading  and  roll  measuring),  with  the  exact  dimensions  given  in  Figure  6. 


FIGURE  5.  GPS  ANTENNA  ARRAY  INSTALLATION 
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FIGURE  6.  GPS  ANTENNA  ARRAY  GEOMETRY 


Table  3  provides  a  list  of  the  systems  relevant  to  this  report,  giving  their  resolutions, 
data  rates  and  their  expected  attitude  accuracies.  Of  course  there  is  at  this  time  no  real 
specification  for  the  accuracy  of  the  3DF  attitude  measurements,  especially  in  a  dynamic 
environment.  The  general  factors  having  an  effect  on  this  accuracy  are  discussed  in  Section  2 
above,  while  the  expected  accuracy  of  this  specific  installation  is  discussed  in  Section  3.2 
below,  where  0  is  heading,  <j>  is  pitch,  vy  is  roll  and  e  =  10'16. 


TABLE  3.  EXPECTED  SYSTEM  PERFORMANCE 


SENSOR 

RESOLUTION 

(degrees) 

0  v 

DATA 

RATE 

(sec.) 

EXPECTED  RMS  ERROR 
(degrees) 

0  <j>,  V|/ 

3DF 

8 

8 

1.0 

=  0.1 

=  0.15 

INS1 

0.005 

0.003 

0.16 

<0.07 

<0.03 

INS2 

0.005 

0.005 

0.51 

0.10 

0.05 

INS1/DGPS 

8 

8 

0.16 

0.01 

0.01 

INS2/DGPS 

8 

8 

0.51 

0.01 

0.01 
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The  resolution,  e,  of  the  binary  3DF  attitude  data  on  output  (which  uses  double 
precision,  or  64  bits)  is  so  much  better  than  the  accuracy  as  to  be  irrelevant.  The  same  is  true 
of  the  output  from  the  INS/GPS  integration  filters. 

The  expected  rms  errors  given  in  Table  3  for  INS2  are  as  specified  for  the  USAF 
Standard  Form,  Fit  and  Function  (F3)  Medium  Accuracy  Inertial  Navigation  Unit  (SNU  84-1), 
which  is  basically  a  one  nautical  mile  per  hour  unit. 

The  approximate  accuracy  performance  figures  given  for  the  INS/DGPS  integrated 
solution  in  Table  3  are  based  on  the  Kalman  filter  covariance  values.  The  validity  of  this 
covariance  information  has  been  verified  to  a  great  extent  by  simply  verifying  the  error  models 
upon  which  it  is  based.  This  was  done  by  a  detailed  examination  of  all  filter  inputs  and  outputs 
(for  example  by  verifying  the  consistency  of  the  residuals  with  the  residual  covariances). 

3.2  Expected  3DF  Performance 

Since  baselines  were  of  the  same  length  in  both  lateral  and  longitudinal  directions,  as 
seen  from  Figure  6,  the  pitch  and  roll  errors  should  be  statistically  of  the  same  magnitude  (by 
symmetry).  However  they  should  be  larger  than  the  heading  errors  by  a  factor  of  about  1 .5  to 
2,  due  to  the  difference  in  vertical  and  horizontal  satellite  geometry,  as  discussed  in  section  2 
above  (see  equation  2).  The  most  relevant  data  from  the  literature  is  from  references  [6c]  and 
[8],  However  the  reference  [6c]  data  was  very  limited  and  apparently  did  not  have  a  good 
reference  system,  whereas  the  reference  [8]  data  was  collected  on  a  truck  driving  down  a 
straight,  smooth  road,  under  close  to  ideal  GPS  reception  conditions.  Since  this  Endeavour 
data  was  collected  24  hours  a  day  for  seven  days,  on  a  gently  pitching  and  rolling  ship,  with 
some  degree  of  superstructure  masking,  we  can  perhaps  expect  more  spurious  data  and  more 
dynamically  induced  errors.  On  the  other  hand  the  antenna  baselines  used  here  were  about  six 
times  longer  than  those  of  reference  [8],  (but  five  times  shorter  than  those  of  reference  [6])  so 
that  on  balance  the  nominal  errors  should  be  roughly  in  line  with  the  corresponding  rows  (6c 
and  8)  of  Table  2:  around  the  0.1  to  0.2  degree  level. 
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FIGURE  7.  VIEW  OF  ANTENNA  ARRAY 


3.3  Reference  System  Performance 

As  mentioned,  van  jus  sensors  were  available  to  provide  reference  data.  The  INSl 
inertial  navigation  system  by  itself  should  provide  sufficiently  accurate  pitch  and  roll  reference 
data,  as  indicated  by  Table  3.  For  heading  reference  however  the  raw  INSl  may  not  be 
adequate.  In  this  case  the  INS1/DGPS  data  may  be  required,  which  is  certainly  accurate 
enough. 

Another  consideration  is  that  the  INSl  (and  INS1/DGPS)  accuracies  listed  in  Table  3 
actually  apply  to  the  discrete  data  points  at  their  time  of  validity.  This  data  was  taken  from  the 
digital  interface  of  INSl  at  a  0.16  second  rate,  which  was  the  slowest  recording  rate  available 
Since  this  reference  data  was  not  synchronized  with  the  one  second  GPS  attitude  times, 
comparison  of  these  two  discrete  data  sets  required  some  data  registration,  to  match 
corresponding  data.  This  data  registration  process  introduces  two  additional  sources  of  error 
latency  error  and  interpolation  error. 

It  could  be  mentioned  that  the  INSl  data  was  available  at  considerably  higher  rates 
from  the  digital  interface  (and  essentially  continuously  from  the  synchro  outputs).  However 
this  was  not  considered  necessary,  as  explained  below,  and  with  limited  data  recording  capacity 
on  a  seven  day  trial,  was  not  practical. 

3.3.1  Latency  Error 

Precise  absolute  time  tagging  of  all  sea  trial  reference  data,  and  in  particular  INSl  and 
INS2  data,  was  accomplished  by  using  GPS  clocks  with  IRIG  B  output  These  clocks  allow 
the  recording  systems  to  time  tag  the  data  with  the  time  of  receipt  of  the  first  byte  of  each 
message  to  an  accuracy  of  several  microseconds.  However  there  was  some  uncertainty  in  the 
exact  time  of  validity  of  the  data  (as  opposed  to  its  time  of  recording).  Even  after  the 
recording  times  are  corrected  for  the  known  latencies  (which  for  the  INSl  reference  data  were 
on  the  order  of  tens  of  milliseconds)  there  was  still  a  time  uncertainty  of  a  few  milliseconds.  (It 
should  perhaps  be  mentioned  that  the  synchro  output  from  INSl,  which  has  negligible  latency, 
could  have  been  used  if  necessary.)  The  attitude  error  introduced  by  this  timing  error  will  be 
determined  by  the  angular  rate  of  the  platform. 

For  the  attitude  reference  system  (INSl)  for  example,  the  uncertainty  in  latency  of  the 
attitude  data  is  St  where: 

St  <  0.002  seconds  (3) 

During  this  time  the  attitude  0  will  have  changed  by  50,  where. 

60  =  6t50/5t  (4) 

where  30/dt  is  the  time  rate  of  change  of  the  attitude  angle.  Thus  this  uncertainty  in  timing  St 
will  cause  an  attitude  error  50  ,  as  given  by  equation  4,  which  is  determined  by  the  platform 
angular  rate  <90/3t.  To  place  a  limit  on  this  latency  error  we  must  place  a  limit  on  this  angular 
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rate. 


As  discussed  in  Appendix  A,  (equation  A8)  the  dynamics  of  the  platform  during  this 
trial  were  such  that  a  maximum  angular  rate  of  less  than  30/dt  =  4.2  degrees/second  is 
expected.  Substituting  in  equation  4  this  indicates  a  data  latency  error  bound  of 

1 60 1  <  0.008  degrees  (5) 

which  is  not  significant. 

Unfortunately  the  3DF  data  itself  was  time  tagged  internally,  and  not  by  the  data 
collection  device  (i.e.  not  at  its  output  times),  so  that  nothing  could  be  said  about  its  data 
latency,  which  is  important  for  use  by  a  dynamic  control  system. 


3.3.2  Interpolation  Error 

Since  the  reference  data,  recorded  at  0.16  second  intervals,  is  not  synchronized  with 
the  one  second  3DF  data,  it  was  necessary  to  interpolate  the  reference  data,  introducing  an 
interpolation  error.  With  linear  interpolation  we  estimate  the  angle  0  to  be 

§(f  +  Af)  =  0(f)  +  Af  x  —(f)  (6) 

dt 

where  in  this  case 


-0.0S  seconds  <  At  <  0.08  seconds  (7) 

The  magnitude  of  the  resulting  attitude  error  j  50 1  will  be  approximately  that  of  the  next  term 
in  the  Taylor  expansion,  namely  At2/2  times  he  angular  acceleration  c^Q/dt2  . 


Af2 

^0 

2 

dt2 

(8) 


Now  the  magnitude  of  the  interpolation  interval  I  At  I  is  at  most  0.08  seconds  in  this 
case.  Thus  to  place  a  limit  on  this  error  it  is  necessary  to  examine  the  angular  acceleration  of 
the  platform.  From  Appendix  A  (equation  A9)  we  see  that  the  magnitude  of  the  angular 
acceleration  is  less  than  4.4  degrees/sec.2.  Thus  from  equation  8  we  see  that  the  interpolation 
error  can  be  bounded  by: 

1 50 !  <  (0.08  sec)2(0.5  )(4.4  degrees/sec2) 

=  0.014  degrees  (9) 

which  is  also  not  very  significant. 
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3.3.3  Total  Reference  System  Error 

When  considering  the  ultimate  accuracy  of  the  reference  data  used,  it  is  important  to 
note  that  the  latency  and  interpolation  errors  are  both  going  to  vary  sinusoidally  in  magnitude, 
and  are  90  degrees  out  of  phase,  since  o*-*  is  proportional  to  angular  rate  and  the  other  to 
angular  acceleration.  Thus  the  bound  on  the  sum  of  these  two  errors  is  much  less  than  the  sum 
of  the  bounds. 

As  already  mentioned,  another  option  if  a  more  accurate  reference  or  more  confidence 
in  the  reference  data  were  required,  would  be  to  run  the  raw  INS1  or  INS2  data  through  a 
Kalman  filter  with  the  differential  GPS  position  measurements.  This  would  improve  the 
attitude  estimates  as  well  as  provide  error  covariance  information  indicating  quantitatively 
what  the  accuracy  is,  as  a  function  of  time.  This  would  also  provide  residual  and  residual 
covariance  information  to  verify  that  the  INS  errors  during  this  sea  trial  did  in  fact  conform  to 
the  error  models  upon  which  the  Kalman  filter  is  based  (and  upon  which  the  INS  performance 
figures  of  Table  3  are  based).  This  would  not  by  itself  remove  the  small  errors  due  to  INS  data 
latency  uncertainty  or  to  interpolation,  and  so  the  performance  would  be  limited  to  about  0.005 
degrees,  as  shown  in  Table  3.  A  higher  order  interpolation  method  however  could  have 
reduced  those  errors,  but  since  they  were  not  significant  this  was  not  done. 

In  fact  the  INSl/DGPS  filter  was  run,  to  verify  that  the  raw  INS1  pitch  and  roll  data 
was  sufficiently  accurate  (which  it  was)  and  to  provide  corrections  to  the  raw  INS1  heading, 
so  that  it  would  also  be  adequate  as  reference  for  the  3DF  data. 

Another  important  aspect  of  the  reference  system  errors  is  their  behaviour  in  time  in 
terms  of  the  associated  autocorrelation  function.  Inertial  system  errors  are  quite  well 
understood  (see  for  example  reference  [13])  and  their  temporal  behaviour  is  distinctly  different 
from  the  behaviour  of  the  errors  attributed  here  to  the  3DF  data.  Because  of  this,  any  inettial 
errors  which  are  present  will  not  contribute  significantly  to  the  error  model  parameters 
attributed  to  the  3DF  measurements.  Neither  will  any  (constant)  data  latency  error,  for  the 
same  reason,  as  will  be  explained  in  Section  6  below. 
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4.  RAW  DATA 


The  3DF  data  set  consisted  of  520,255  data  points,  collected  at  a  one  Hz  rate  over  the 
duration  of  the  trial,  which  was  586,800  seconds  long  (almost  7  days).  Thus  there  was  1 1%  of 
the  trial  for  which  there  was  no  3DF  data,  primarily  because  of  gaps  in  the  GPS  coverage 
(agreeing  closely  with  the  1 1%  missing  differential  GPS  reference  data).  There  were  also  brief 
periods,  once  a  day,  when  the  3DF  data  recording  was  shut  down  to  close  off  data  files 

A  further  10%  of  the  available  data  was  not  "valid"  Figure  8  shows  the  attitude  reset 
flag  values  throughout  the  trial,  as  received  from  the  3DF  receiver.  This  flag  is  "0"  if  the  phase 
ambiguities  have  been  resolved  and  a  carrier  phase  attitude  has  been  computed.  It  is  "1"  if  the 
phase  ambiguity  has  not  been  solved,  in  which  case  the  pitch  and  roll  are  set  to  zero  and  a 
codephase  (as  opposed  to  carrier  phase)  heading  estimate  is  given.  For  purposes  of  this  report, 
only  carrier  pha*a  attitude  data  has  been  analyzed.  This  "valid”  data  constituted  about  90%  of 
the  available  data  (469,590  of  the  one  second  3DF  data  points  were  flagged  ”0",  out  of  the 
total  520,255  3DF  data  points).  Thus  the  breakdown  is  as  follows: 

Total  trial  period  586,800  seconds  (100%) 

Total  3DF  data  520,255  seconds  (89%) 

Valid  3DF  data  469,590  seconds  (80%) 

Figures  9  to  1 1  show  the  raw  inertial  heading,  pitch  and  roll  respectively,  throughout 
the  trial.  These  plots  are  included  simply  to  provide  an  indication  of  the  dynamics  experienced 
during  this  trial  The  trajectory  plot  of  Figure  4  also  gives  an  indication  of  the  various 
headings  assumed. 

Figures  12  to  14  show  the  3DF  attitude  (at  1  Hz)  in  comparison  to  the  reference 
attitude  (at  6.25  Hz),  for  a  brief  (100  second)  interval.  These  illustrate  the  short  term 
dynamics,  due  largely  to  wave  motion,  and  at  the  same  time  give  a  first  impression  of  the 
quality  of  the  3DF  data.  An  appreciation  of  these  short  term  dynamics  is  important  when 
analyzing  and  interpreting  the  performance  of  the  3DF  data,  and  these  plots  will  be  referred  to 
again  below 

Figure  15  shows  the  MRMS,  which  is  the  average  double  difference  carrier  phase 
residual,  as  supplied  by  the  3DF  receiver.  This  is  supposed  to  be  typically  2  to  3  millimeters 
Thus  Figure  15  seems  to  indicate  that  there  were  no  unusual  problems  with  these  phase 
measurements. 

Figure  16  shows  the  BRMS,  which  is  the  RMS  error  for  the  differences  between 
calibrated  baseline  magnitudes  and  computed  baseline  magnitudes  for  the  three  vectors  formed 
by  Antenna  l  to  the  other  three  antennas.  This  is  supposed  to  be  typically  1  to  3  centimeters, 
increasing  under  high  PDOP  conditions.  Thus  Figure  16  seems  to  indicate  that  there  were  no 
serious  problems  with  PDOP  or  with  these  calibration  values. 
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FIGURE  8a:  AVAILABILITY 
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FIGURE  15:  3DF  MRMS 


5.  MEASURED  3DF  ERROR  BEHAVIOUR 


During  this  trial  there  were  14  Block  II  GPS  satellites  plus  4  Block  I  satellites 
operational.  The  full  GPS  operational  capability  is  expected  in  1995,  at  which  time  the  GPS 
constellation  will  consist  of  24  (Block  II)  satellites.  Since  the  constellation  was  not  fully 
operational  at  the  time  of  this  trial,  some  brief  periods  of  poor  satellite  geometry  were 
experienced,  causing  some  abnormally  large  errors.  Indeed  about  9.7%  of  the  3DF  attitude 
data  was  flagged  as  invalid,  as  illustrated  in  Figure  8.  This  data  was  therefore  removed  and  all 
statistical  results  shown  here  are  based  on  the  remaining  data. 

Since  this  report  is  dealing  with  the  real-time  dynamic  performance  of  a  system  which 
only  provides  discrete  output,  care  must  be  taken  to  define  what  errors  are  being  examined. 
For  example  data  interpolation  is  not  possible  in  real  time.  Therefore  to  indicate  the  potential 
3DF  performance,  we  have  attempted  to  separate  the  actual  measurement  errors  from  the 
errors  due  to  changes  in  the  attitude  between  measurements. 

We  therefore  first  look  only  at  the  errors  in  the  3DF  measurements  at  the  one  second 
marks.  These  one  second  errors  will  be  referred  to  as  the  discrete  dynamic  errors.  These 
can  be  interpreted  as  the  errors  seen  by  a  user  who  did  not  require  continuous  attitude  data, 
such  as  an  integrated  system  with  an  inertial  component.  They  could  also  be  seen  as  an  upper 
bound  on  the  potential  performance  of  such  an  integrated  system  which  used  a  very  low  cost 
INS  (or  AHRS).  However  these  figures  can  be  very  misleading  for  real-time  dynamic 
applications,  as  discussed  below. 

These  discrete  dynamic  errors  however  do  not  account  for  the  very  real  dynamic 
limitations  in  the  one  second  measurements.  For  example,  if  attitude  data  was  being  used  in 
real-time  at  a  50  Hz  rate  but  was  only  updated  once  a  second,  then  any  change  in  attitude  over 
the  one  second  update  intervals  (due  to  platform  dynamics)  would  effectively  be  a  real  error  in 
the  estimate,  even  it  were  perfectly  accurate  at  the  one  second  marks.  This  effect  can  be  seen 
Figures  12  to  14,  in  the  case  of  the  3DF  system  on  a  ship  at  sea.  Here  the  platform  dynamics 
between  3DF  samples  can  be  clearly  seen.  Since  the  reference  data  was  recorded  at  a  higher 
rate  (about  6  Hz),  a  more  honest  evaluation  of  the  continuous  3DF  performance  is  found  by 
comparing  the  single  3DF  value  for  each  one  second  interval  with  the  reference  values  over  the 
entire  interval.  Errors  found  in  this  way  will  be  referred  to  as  the  continuous  dynamic 
errors.  Since  this  accurately  reflects  the  real-time  dynamic  performance,  it  depends  strongly  on 
the  platform  dynamics  and  the  measurement  update  rate.  The  results  found  here  therefore 
apply  only  to  a  particular  platform  (small  ship)  and  rate  (1  Hz),  and  are  intended  primarily  to 
illustrate  how  misleading  the  discrete  dynamic  performance  can  be. 

5.1  Discrete-Dynamic  Performance 

The  discrete  dynamic  errors  of  the  3DF  heading,  pitch  and  roll  are  shown  in  Figures  17 
to  19,  where  the  time  histories  throughout  the  trial  are  shown.  Only  valid  carrier  phase  derived 
attitude  data  is  shown.  From  these  plots  it  is  evident  that  the  pitch  and  roll  errors  are 
dominated  by  constant  biases,  presumably  due  to  misalignment. 
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Figures  20  to  22  show  the  same  data  on  a  smaller  scale,  so  that  the  nominal  data  can  be 
more  easily  seen.  From  these  it  can  be  seen  that  the  pitch  and  roll  errors  are  dominated  by 
noise,  while  the  heading  error  has  a  more  complex  behaviour,  with  some  drifting  components 
and  sudden  shifts  in  "bias".  Appendix  B  presents  this  data  on  an  expanded  time  scale  for 
clearer  viewing,  showing  daily  plots.  This  appendix  also  shows  the  heading  errors  with  some 
smoothing,  so  that  the  drifting  errors  and  changes  in  the  "bias"  are  not  obscured  by  the  noise. 

The  statistics  derived  from  this  data  are  presented  in  Table  4.  Since  the  pitch  and  roll 
errors  are  obviously  dominated  by  biases,  presumably  due  to  installation/calibration  errors, 
statistics  have  also  been  derived  from  data  with  these  biases  removed,  to  better  illustrate  the 
stability  and  the  potential  accuracy  performance. 

As  a  cautionary  note,  it  should  be  mentioned  that,  since  the  errors  obviously  have  a 
large  component  of  noise,  the  amount  of  data  compressed  onto  the  plots  of  Figures  1 7  to  22 
makes  it  difficult  to  visually  judge  the  rms  values  or  standard  deviations.  With  a  plotting 
resolution  of  300  dots  per  inch,  presenting  469,577  data  point  against  an  8  inch  scale  means 
that  at  each  time  pixel  about  200  data  points  are  superimposed,  so  that  only  the  two  extreme 
points  from  each  set  are  visible.  This  of  course  has  a  magnifying  effect  on  the  apparent  size  of 
the  errors. 

Figures  23  to  25  present  a  detailed  view  of  the  discrete  attitude  errors,  for  the  time 
period  corresponding  to  Figures  12  to  14. 


TABLE  4.  DISCRETE-DYNAMIC  3DF  ATTITUDE  ERRORS 


RMS 

STANDARD 

PERCENTILE  (degrees) 

ERROR 

DEVIATION 

raw 

bias  removed 

SOURCE 

(degrees) 

(degrees) 

68% 

95% 

68% 

95% 

100% 

HEADING 

0.47 

0.47 

0.08 

0.13 

0.04 

0.08 

56.7 

PITCH 

1.14 

0.27 

1  13 

1.20 

0.05 

0.11 

4.5 

ROLL 

1.28 

0.26 

1.27 

1.34 

0.05 

0.11 

6.5 

The  rms  and  standard  deviation  numbers  in  Table  4  are  based  on  469,577  samples, 
taken  at  a  one  Hz  rate  over  seven  days  (all  available  carrier  phase  based  attitude 
measurements)  The  percentiles  are  based  on  samples  taken  every  20  seconds  over  the  seven 
days. 

Figures  26  to  28  show  the  complete  set  of  percentile  values  for  the  valid  data  (with 
biases  removed).  This  is  simply  a  plot  of  all  error  magnitudes,  sorted  by  magnitude.  It  can  be 
interpreted  as  a  sample  cumulative  probability  density  function  as  follows.  If  the  values  of  the 
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(10) 


percentile  function  for  1 50 1  is  V(p)  for  percentiles  0<p<  100,  then 
Prob{  1 60 1  <;  V}  =  p/100 

Also  noteworthy  is  the  fact  that  the  rms  and  standard  deviation  statistics  in  Table  4  are 
very  large  when  compared  to  the  percentile  values,  for  normally  distributed  data  the  rms  value 
should  correspond  to  the  68  percentile.  As  can  be  seen  from  the  percentile  plots,  the  large  rms 
values  are  due  to  a  very  small  number  of  relatively  large  errors.  If  this  were  filtered  with 
inertial  data  then  these  spurious  points  could  be  eliminated.  In  fact  when  the  worst  1%  of  the 
errors  are  removed,  then  the  standard  deviation  of  the  remaining  data  is  0.045  degrees  in 
heading  and  0.054  degrees  in  pitch  and  roll,  which  very  closely  matches  the  68  percentiles 
(bias  removed),  as  would  be  expected.  Figures  29  to  3 1  show  the  error  percentiles  rescaled,  so 
that  the  more  meaningful  values  can  be  clearly  seen. 

The  percentile  statistics  are  therefore  perhaps  the  most  significant.  These  illustrate  the 
nominal  performance  of  about  0.05  degrees  (or  about  one  milliradian)  at  the  one  sigma  level 
and  about  0. 1  degree  (or  about  two  milliradians)  at  the  two  sigma  (95%)  level.  These  discrete- 
dynamic  errors  are  examined  in  more  detail  in  Section  6  below. 


5.2  Continuous-Dynamic  Performance 

For  real-time  use,  as  a  stand-alone  system  on  a  moving  platform,  the  continuous- 
dynamic  performance  is  probably  the  most  relevant.  The  continuous  dynamic  errors  of  the 
3DF  heading,  pitch  and  roll  are  shown  in  Figures  32  to  34,  where  the  time  histories  throughout 
the  trial  are  shown.  These  results  are  of  course  highly  dependent  on  the  platform  dynamics, 
which  in  this  case  were  quite  benign,  and  the  measurement  rate,  which  at  1  Hz  was  quite  slow 
As  seen  in  Figures  9  to  14,  the  ship  experienced  fairly  low  amplitude  (less  than  about  4 
degrees)  attitude  oscillations  with  about  a  5  to  10  second  period. 

As  with  the  discrete  data,  only  valid  carrier  phase  derived  attitude  data  is  shown. 
Figures  35  to  37  show  the  same  data  on  a  smaller  scale,  so  that  the  nominal  data  can  be  more 
easily  seen.  The  statistics  derived  from  this  data  are  presented  in  Table  5.  Since  the  pitch  and 
roll  errors  are  obviously  dominated  by  biases,  presumably  due  to  installation/calibration  errors, 
statistics  have  also  been  derived  from  data  with  these  biases  removed.  Comparing  these  results 
with  the  discrete-dynamic  statistics  of  Table  4,  it  becomes  clear  that  the  continuous-dynamic 
errors  are  dominated  by  the  "extrapolation  error"  between  the  one  second  measurements. 

Figures  38  to  40  present  a  detailed  view  of  these  attitude  errors,  for  the  time  period 
corresponding  to  Figures  12  to  14.  These  illustrate  the  significant  error  growth  between  one 
second  measurement  updates,  due  to  platform  dynamics.  Clearly  an  increased  update  rate,  if 
possible,  and/or  a  careful  extrapolation  method,  would  be  desirable. 

Figures  41  to  43  show  the  complete  set  of  percentile  values  for  the  valid  data  (with 
biases  removed).  As  with  the  discrete  data,  the  statistics  here  are  skewed  by  a  very  small 
number  of  bad  data  points  (less  than  1%).  Thus,  as  with  the  discrete  results,  the  percentiles  are 
perhaps  the  most  significant  statistics  to  quote  and  Figures  44  to  46  show  these  error 
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percentiles  rescaled,  so  that  the  more  meaningful  values  can  be  clearly  seen. 


TABLE  5.  CONTINUOUS-DYNAMIC  3DF  ATTITUDE  ERRORS 


RMS 

STANDARD 

PERCENTILE  (degrees) 

ERROR 

DEVIATION 

raw 

bias  removed 

SOURCE 

(degrees) 

(degrees) 

68% 

95% 

68% 

95% 

100% 

HEADING 

0.62 

0.61 

0.13 

0.37 

0.13 

56.1 

PITCH 

1.24 

0.56 

1.21 

1.90 

0.28 

0.91 

8.1 

ROLL 

1.35 

0.50 

1.34 

1.93 

0.26 

0.84 

7.8 

The  rms  and  standard  deviation  numbers  in  Table  5  are  based  on  2,974,775  samples, 
taken  at  0. 16  second  intervals  over  seven  days.  As  with  the  discrete  results,  the  percentiles  are 
based  on  every  20th  sample  taken,  (which  in  this  case  is  every  3.2  seconds)  over  the  seven 
days. 

Since  the  characteristics  of  the  continuous  errors  can  be  derived  from  those  of  the 
discrete  errors  if  the  platform  dynamics  are  known,  these  continuous  errors  will  not  be  further 
examined. 
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FIGURE  33:  COKT1NOUS  PITCH  ERROR  VS  TIME 
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FIGURE  37:  CONtlNOUS  ROLL  ERROR  VS  TIME 
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FIGURE  42:  CONTINUOUS  PITCH  ERROR  PERCENTILES 


6.  STOCHASTIC  ERROR  MODEL 


In  this  section  a  specific  stochastic  error  model  will  be  developed  for  each  of  the  3DF 
attitude  error  components.  Besides  offering  insight  into  the  error  behaviour  of  the  3DF 
attitude  measurements,  these  models  will  be  in  a  form  suitable  for  use  by  a  Kalman  filter.  As 
expected,  these  errors  will  be  shown  to  have  characteristics  which  are  complementary  to  the 
error  characteristics  of  an  inertial  system,  making  a  Kalman  filter  integration  of  the  two  an 
attractive  option. 

The  stochastic  models  normally  used  by  Kalman  filter  designers  to  describe  random 
variables  x(t),  are  relatively  simple  linear  models  such  as  random  biases,  white  noise  processes, 
first  order  Markov  processes,  random  walks  or  periodic  processes.  These  are  adequately 
described  in  reference  [14],  They  are  usually  distinguished  by  their  distinctive  autocorrelation 
functions  cp(A/)-  The  autocorrelation  function  cp^A/),  of  a  random  variable  x(t),  is  defined  to 
be  the  expected  correlation  between  values  of  x(/)  separated  in  time  by  At  seconds: 

cPxx(A/)  ■  E{x(t)x(t+At)}  (11) 

Perhaps  the  most  useful  stochastic  model  is  the  first  order  Markov  process,  since  the 
random  bias  and  white  noise  are  both,  in  some  sense,  a  special  case  of  this  Markov  process. 
The  Markov  process  x(t)  is  described  by  two  parameters:  its  root  mean  square  value  a  and  its 
autocorrelation  time  x.  It  has  the  following  autocorrelation  function: 

<Pxx(A/)  =  a2e*A//t  (12) 

For  obvious  reasons  this  type  of  process  is  often  referred  to  as  an  "exponentially 
correlated  process".  If  the  correlation  time  x  is  made  very  large,  this  Markov  process  will  be 
so  strongly  correlated  that  it  will  essentially  behave  as  a  random  bias,  and  if  the  correlation 
time  x  is  made  very  small  (much  smaller  than  the  sample  time),  then  it  will  be  so  weakly 
correlated  that  it  will  essentially  behave  as  uncorrelated  white  noise.  Equation  12  therefore 
provides  an  ideal  "template"  to  use  in  extracting  the  error  model  parameters  (a  and  x)  from  an 
autocorrelation  function. 

Now  the  definition  given  by  equation  11  can  be  used  to  obtain  a  "sample 
autocorrelation  function"  from  the  measured  x(t)  (the  3DF  attitude  errors  in  this  case)  Then 
by  matching  the  plot  of  this  autocorrelation  function  (p**  to  the  template  of  equation  12,  the 
model  parameters  o  and  x  can  be  easily  extracted  as  follows.  The  initial  value  is  o: 

a2  -  tp^O)  (13) 

and  the  point  at  which  (p^  drops  to  1/e  of  its  initial  value  is  x: 

<MT)  =  (14) 

or  alternately,  the  slope  of  cp^  at  the  origin  is  -cp^OJ/x: 
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Therefore,  in  order  to  match  the  3DF  errors  to  the  appropriate  model(s),  it  is  useful  to 
examine  the  sample  autocorrelation  functions. 


6.1  Autocorrelation  Functions 

As  described  in  reference  [14],  the  sample  autocorrelation  function  (p**  of  the  variable 
xn,  n=l...N  (discrete  samples  of  a  stochastic  process  x(t)  at  times  nAi  where  x(t)  may  be  a 
vector)  is  defined  to  be: 


N-n 


|  —n 

=  — - -  £(*,  -«)(*,>„  ~m) 

N -tt- 1 


/t=0,1,...,JV-2 


(16) 


i=l 


where  m  is  the  sample  mean: 


m 


1  N 

=  i-Yx, 


(17) 


where  it  can  be  seen  from  equation  16  that,  as  expected,  tp^fO)  is  the  mean  square  value.  The 
relationship  between  this  sample  autocorrelation  function  and  the  true  autocorrelation  function 
is  discussed  in  reference  [14]  (page  89),  where  it  is  indicated  that  the  sample  function  is  a  very 
good  approximation  of  the  true  function  provided  that  the  sample  length  is  much  longer  than 
the  correlation  time  of  the  process.  In  fact,  for  a  Markov  process,  the  uncertainty  in  the  mean 
o  and  correlation  time  x  ,  as  derived  from  the  sample  autocorrelation  function  is: 


For  small  At/x  («  1),  as  will  usually  be  the  case,  this  can  be  simplified  to: 


(20) 
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In  this  case  the  sample  period  is  At  =  1  and  N  =  469,577.  Therefore  if  t  s  600 
seconds  and  o  =  .04  degrees  (and  if  the  process  really  is  a  first  order  Markov  process)  then 
using  the  sample  autocorrelation  function  to  estimate  a  and  x  introduces  only  about  0  00015 
degrees  (4%)  error  in  the  case  of  o  and  30  seconds  (5%)  in  the  case  of  x ,  which  are  not 
significant  for  most  purposes. 

For  modeling  purposes  it  is  most  appropriate  to  first  model  the  discrete-dynamic  errors, 
since  the  continuous-dynamic  error  model  can  be  obtained  from  this  by  simply  adding  the 
expected  errors  due  to  platform  dynamics.  The  complete  sample  autocorrelation  functions  of 
the  discrete  attitude  errors,  (with  the  biases  removed  and  the  worst  1%  of  errors  removed)  are 
therefore  shown  in  Figures  47  to  49. 

Now  any  values  of  (p  less  than  the  square  of  the  rms  accuracy  of  the  reference  system 
are  not  significant  and  should  be  ignored.  In  this  case,  based  on  INS/DGPS  reference  accuracy 
from  Table  3,  any  correlation  values  below  about  0.0001  degrees2  should  be  ignored.  Based 
on  INS1  reference  data  alone,  these  limits  arc  much  higher,  and  recourse  must  be  taken  to  the 
known  correlation  characteristics  of  the  INS.  However,  from  Figures  47  to  49  it  can  be  seen 
that  there  is  no  significant  heading  error  autocorrelation  for  nAt  >  10,000  seconds  and  no 
significant  pitch  and  roll  error  autocorrelation  for  nAt  >  3,000  seconds. 

Figures  50  to  52  show  the  first  8,000  seconds  of  these  functions,  from  which  it  can  be 
seen  that  there  is  at  least  some  temporally  correlated  error,  and  in  the  case  of  heading  a 
substantial  amount.  In  fact,  ignoring  for  the  moment  the  smaller  structure  beyond  8,000 
seconds  (which  will  be  briefly  discussed  below),  these  sample  autocorrelation  function  plots 
appear  to  be  very  closely  matched  to  the  correlation  function  of  a  first  order  Markov  process 
(exponentially  decaying  as  in  equation  12).  Closer  examination  however,  reveals  that  there  are 
strong  peaks  at  At  =  0,  as  is  perhaps  more  apparent  in  Figures  53  to  55  which  show  the  first 
800  seconds.  The  small  fine  structure  at  0  <  At  <  50  seconds  will  also  be  explained  below 
For  the  moment  we  will  focus  on  modeling  the  dominant  error  components,  the  white  noise 
and  the  Markov  processes. 

If  each  of  these  3DF  attitude  errors  (heading,  pitch  and  roll)  is  assumed  to  be  the  sum 
of  two  stationary  processes: 

56  =  56,  +  S62  (21) 

where  56  [  is  white  noise  and  S62  is  a  Markov  process,  then  if  these  are  independent 
( uncorrelated  with  each  other)  it  can  be  shown  that  the  autocorrelation  function  for  56  is 
simply  the  sum  of  the  autocorrelation  functions  of  50 1  and  562.  Of  course  this  generalizes  to 
the  sum  of  n  independent  processes,  so  that  the  autocorrelation  function  can  be  used  to 
decompose  the  error  into  component  parts.  It  should  be  noted  that  if  errors  come  from 
independent  sources,  such  as  for  example  receiver  noise  and  multipath,  then  they  are  almost 
certainly  independent. 

Since  the  autocorrelation  function  of  a  white  noise  process  is  zero  except  at  At  =  0,  the 
peaks  at  At  =  0  are  due  to  the  white  noise  components  and  can  be  removed  to  find  the  model 
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parameters  for  the  Markov  component.  The  magnitude  of  the  peaks  above  the  exponentially 
decaying  portion,  in  Figures  47  to  55,  are  therefore  the  mean  square  of  the  white  noise 
component.  Unfortunately  these  peak  values  at  At  =  0  cannot  actually  be  seen  from  the  data 
plots  (since  they  drop  off  so  sharply  at  At  =  1,  which  at  these  scales  is  not  clearly  visible). 
However  the  numerical  values  can  of  course  be  easily  obtained  from  the  corresponding  data 
files.  This  error  model  parameter  extraction  process  is  discussed  below,  and  the  results  are 
summarized  in  Table  6. 

Pitch 

From  the  sample  autocorrelation  data  file,  the  total  a2  for  the  pitch  error  is  cp^(O)  = 
0.0028  degrees2.  From  Figure  54  the  Markov  a22  can  be  estimated  at  0 .0012  degrees^.  This 
puts  the  white  noise  cr^2  estimate  at  0.0016  degrees2.  The  Markov  correlation  time  is  the 
point  where  the  Markov  <p  drops  below  <p^(0)/e  =  0.0012/2.718  =  0.00044,  which  appears  to 
be  at  about  t2  =  620  seconds. 

Roll 

From  the  sample  autocorrelation  data,  the  total  a2  for  the  roll  error  is  <pw(0)  =  0  0028 
degrees2.  From  Figure  55  the  Markov  o22  can  be  estimated  at  0.0012  degrees2  which  puts  the 
white  noise  Oj2  estimate  at  0.0016  degrees2.  The  Markov  correlation  time  is  the  point  where 
the  Markov  <p  drops  below  <pVVJ/(0)/e  =  0.0012/2.718  =  0.00044,  which  appears  to  be  at  about 
x2  =  550  seconds. 

Heading 

The  heading  error,  as  seen  in  Figure  47,  has  a  somewhat  more  complicated 
autocorrelation  function,  as  should  perhaps  be  expected  (and  as  will  be  discussed  below).  In 
this  case  there  seems  to  be  a  dominant  longer  correlation  time  Markov  process  plus  several 
smaller  processes,  none  of  which  can  be  considered  significant.  (They  consist  of  small  periodic 
processes  with  12  hour,  24  hour  and  84  minute  periods,  and  perhaps  a  very  small,  short 
correlation  time  Markov  process  (these  last  two  are  only  visible  by  close  examination  of  Figure 
50).  Because  these  are  of  questionable  statistical  significance,  only  the  white  noise  and  long 
term  Markov  process  will  be  modeled  here. 

From  the  sample  autocorrelation  data,  the  total  o2  for  the  heading  error  is  <Peo(0)  = 
0.0013  degrees2.  From  Figure  53  the  Markov  o22  can  be  estimated  at  0.0005  degrees2  which 
puts  the  white  noise  a j2  estimate  at  0.0008  degrees2.  The  Markov  correlation  time  is  the 
point  where  the  Markov  <p  drops  below  cpee(0)/e  =  0.0005/2.718  =  0.0002,  which  appears  to 
be  at  about  t2  =  4,000  seconds.  To  see  this,  one  should  remove  the  very  small  "Schuler  bump" 
(due  to  inertial  errors)  at  about  84  minutes  (5,040  seconds)  by  matching  a  straight  line  to  the 
sample  autocorrelation  function  on  both  sides  of  this  bump.  According  to  equation  19,  there 
should  be  an  uncertainty  in  t2  of  about  500  seconds. 

The  biases,  which  have  been  removed  from  the  data  prior  to  these  calculations,  are 
installation  dependent  and  therefore  more  difficult  to  model  in  this  way.  Since  only  one 
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installation  was  done  for  this  trial,  the  bias  model  parameters  in  Table  6  are  based  on  a  sample 
of  one  and  represent  an  order-of-magnitude  estimate  only.  They  must  therefore  be  treated 
carefully. 

It  should  also  be  kept  in  mind  that  these  models  are  based  on  data  from  which  the 
spurious  data  (worst  1%  of  errors)  have  been  removed.  While  this  spurious  data  had  a 
significant  effect  on  the  performance  statistics  (as  shown  in  Table  4  above)  they  could 
presumably  be  rejected  by  an  integrated  system.  In  a  stand  alone  GPS  system  however,  these 
spurious  points  cannot  be  ignored.  It  should  also  be  kept  in  mind  that  these  discrete  errors  do 
not  include  the  latency  or  interpolation  errors  which  would  be  present  (and  very  significant)  in 
the  continuous  dynamic  situation  with  a  stand  alone  GPS  system. 


TABLE  6.  DISCRETE  3DF  ERROR  MODEL  PARAMETERS 


BIAS* 

WHITE 

NOISE 

MARKOV 

<*0 

(degrees) 

(degrees) 

o2 

(degrees) 

*2 

(seconds) 

HEADING 

50 

0.1 

0.03 

0.022 

4,000 

PITCH 

54 

1.0 

0.04 

0.035 

600 

ROLL 

5vy 

1.0 

0.04 

0.035 

600 

*  order  of  magnitude  only 


A  simple  interpretation  of  this  model  would  be  that  the  Markov  components  were  due 
to  multipath  errors,  the  white  noise  components  due  to  receiver  noise  and  the  bias  due  to 
installation  alignment  (calibration)  error. 

The  receiver  noise,  which  should  introduce  the  same  level  of  phase  error  in  each 
measurement,  should  therefore  produce  a  slightly  smaller  heading  error,  due  to  the  geometric 
effect  mentioned  in  Section  2.  This  is  consistent  with  the  measured  noise  components  shown 
in  Table  6. 

The  multipath  error,  on  the  other  hand,  may  not  introduce  the  same  phase  error  in  each 
measurement.  In  fact  it  seems  plausible  that  low  satellite  elevation  angle  signals,  which 
contribute  primarily  to  the  heading  measurement,  would  be  much  more  susceptible  to 
multipath  errors  than  high  angle  signals,  which  contribute  primarily  to  the  pitch  and  roll.  This 
is  consistent  with  the  measured  Markov  components  shown  in  Table  6.  There  in  fact  seems  to 
be  strong  circumstantial  evidence  that  a  large  portion  of  the  correlated  heading  error  is  due  to 
multipath.  There  are  various  times  where  a  sudden  shift  in  heading  error  takes  place, 
corresponding  exactly  with  a  sudden  change  in  ship  heading,  as  seen  at  time  40,000  sec.  in 
Figure  Bl,  time  26,000  sec.  in  Figure  B2,  time  70,000  sec.  in  Figure  B5  and  time  46,000  sec. 
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in  Figure  B6.  Furthermore,  the  4,000  second  correlation  time  of  the  heading  Markov  error 
corresponds  roughly  to  the  expected  correlation  time  of  the  multipath  error,  which  changes 
with  satellite  geometry.  The  satellites,  in  12  hour  orbits,  move  through  about  33  degrees  in 
4,000  seconds,  which  corresponds  roughly  to  a  "significant"  change  in  geometry  (enough  to 
substantially  change  the  multipath  errors). 

Because  of  these  maneuver-induced  discontinuities  in  the  heading  error,  this  multipath 
error  does  not  truly  behave  as  a  pure  Markov  process.  These  discontinuities  could  be  modeled 
as  a  Poisson  process  (see  reference  [15])  which  can  also  be  modeled  within  a  Kalman  filter  as 
described  in  reference  [16],  but  this  degree  of  modeling  fidelity  is  likely  unnecessary.  One 
might  expect  this  error  to  have  a  periodic  component  as  well  as  the  Markov  component,  since 
the  GPS  geometry  repeats  every  12  hours  (provided  the  platform  does  not  move  a  great 
distance).  The  platform  orientation  relative  to  the  GPS  constellation  would  not  generally 
repeat  in  the  same  way,  but  since  this  trial  used  primarily  cardinal  headings  (north,  south,  east 
and  west,  with  some  north-west  and  south-east),  there  was  a  significant  probability  of 
coincident  geometry  from  one  period  to  the  next. 

Another  factor  which  may  affect  heading  differently  than  pitch  and  roll  is  ionospheric 
and  atmospheric  refraction,  which  will  be  much  greater  for  low  elevation  satellites  than  high. 
While  the  effect  on  time  of  arrival  of  the  signal  is  removed  by  the  double  differencing  technique 
used  to  obtain  attitude,  the  (possible)  effect  of  this  refraction  on  the  direction  of  incidence  will 
be  relevant.  The  correlation  function  for  this  type  of  error  would  be  expected  to  also  have  a 
12-hour  periodic  component  plus  a  Markov  component.  A  full  analysis  of  this  effect  is  beyond 
the  scope  of  this  report. 

Various  other  error  sources  could  also  contribute  to  a  lesser  extent,  and  could  be 
modeled  for  some  special  applications.  For  example,  a  data  registration  error  (due  to  an 
unknown  data  latency,  as  discussed  in  Section  3.3  above)  would  introduce  a  periodic  error, 
with  the  period  being  that  of  the  platform  dynamics,  which  was  about  10  seconds  in  this  case. 
This  would  introduce  a  decaying  periodic  component  to  the  autocorrelation  function,  with  the 
same  period.  Indeed  such  a  component  can  be  clearly  seen  in  all  three  of  the  sample 
autocorrelation  functions  in  Figures  53  to  55.  The  small  magnitude  of  these  periodic 
components  (about  0.01  degrees),  when  compared  to  the  observed  3DF  error  components, 
shows  that  the  data  registration  was  adequate  for  this  purpose. 

Another  factor  to  consider  is  the  internal  Kalman  filtering  within  the  3DF  receiver,  as 
mentioned  in  reference  [6].  This  could  have  the  effect  of  reducing  the  white  noise  component 
of  the  error  at  the  expense  of  the  Markov  component.  The  filter  gains  used  will  determine  the 
correlation  time  of  the  resulting  "error"  (higher  gain  leading  to  shorter  correlation  time). 
Hopefully  this  correlation  time  is  very  short. 

Finally  the  errors  in  the  reference  system  must  also  be  carefully  considered.  Reference 
[13]  provides  an  excellent  description  of  inertial  system  error  behaviour.  The  inertial  error 
would  explain  the  24-hour  (86,400  second)  periodic  component  in  the  sample  heading 
autocorrelation  function  as  well  as  the  very  small  Shuler  periodic  components  (5040  seconds) 
in  all  three  sample  autocorrelation  functions. 
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6.2  Cross-Correlation  Functions 


For  use  as  measurements  in  a  Kalman  filter,  the  cross-correlation  is  also  important. 
The  cross-correlation  function  of  two  processes  is  defined  by  the  obvious  generalization  of 
equation  (1 1): 


<Pxy(At)  ^  E{x(/)y(f+A/)}  (22) 

and  the  sample  cross-correlation  function  is  the  corresponding  generalization  of  equation  (16): 


Vxyinto)  3— 7  £(*/ -mx){yi+n~my)1  «= <U,~  JV-2 
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The  Markov  processes  would  be  expected  to  be  cross-correlated  since  they  arise  from 
different  projections  of  the  same  carrier  phase  errors.  (This  will  be  properly  accounted  for  if 
the  phase  errors  are  modeled.)  It  is  the  cross-correlation  of  the  measurement  noise  that  is 
primarily  of  interest,  since  negligible  cross-correlation  here  will  simplify  the  Kalman  filter. 

Figures  56  to  61  show  these  three  cross-correlation  functions  at  different  scales.  From 
these  it  can  be  seen  that  the  three  measurement  noise  components  do  not  have  any  significant 
cross-correlation  beyond  a  few  seconds,  as  would  be  expected.  However,  the  important  error 
model  parameters  here  (from  a  Kalman  filter  perspective)  are  the  cross-covariance  values: 


(Pxy(0)  h  EWOMO} 


(24) 


and  the  correlation  coefficients. 

E{x(/)K/)}-E{x(/)}E(M/)) 


Pxy 


(25) 


From  the  data  shown  in  Figures  59  to  61,  and  using  equation  24  we  have  the  sample  cross- 
covariance  values: 


E{69  64>}  =  0.000023  degrees2 

E{S0  6y}  =  -0.00018  degrees2 

E{  64>  5vj/}  =  -0.00058  degrees2 

which  imply  the  following  correlation  coefficients: 

E{60  54>)  -  E{50}E{64>> 


(26) 


0.000023-0 


Pey 


p<j>y 


VO.  0013x0. 0028 
0.012 

E{SeSy}-E{Se}E{SYt;} 

-0.00018-0 
VO.  0013  x  0.0028 


-0.094 

E{S<t>  8m/}  -  E{5<}>}E{5v|/} 

-0.00058-0 

0.0028 


(27) 


=  -0.21 


These  correlation  coefficients  (which  can  in  general  range  from  -1  to  1)  suggest  that 
there  is  a  significant  degree  of  cross-correlation,  especially  between  the  pitch  and  roll  errors. 


6.3  Kalman  Filter  Model 

With  this  understanding  of  the  error  sources,  one  can  begin  to  formulate  an  error  model 
suitable  for  use  in  a  Kalman  filter  based  integrated  system.  The  biases  can  be  modeled  as 
independent  random  constants  (or  very  long  correlation  time  Markov  processes)  in  the  state 
vector.  The  white  noise  component  is  more  appropriately  modeled  as  measurement  noise, 
outside  the  state  vector.  Whether  the  Markov  components  should  be  modeled  explicitly  in  the 
state  vector  or  more  simply  absorbed  into  the  noise,  is  dependent  on  the  application  and  is 
beyond  the  scope  of  this  report. 

If  these  Markov  states  were  used,  it  may  be  preferable  to  model  them  at  their  source,  as 
independent  carrier  phase  errors,  rather  than  as  attitude  errors,  since  the  carrier  phase  errors 
will  be  independent  to  a  much  greater  extent.  This  is  because  the  phase  errors  from  each 
receiver  channel  (noise  and  multipath  etc.)  will  generally  contribute,  in  vaiying  degrees,  to  all 
three  components  of  the  attitude  error.  This  will  introduce  cross-correlations  into  the  attitude 
errors  which  can  be  properly  accounted  for  only  by  modeling  the  phase  errors.  If  the  GPS 
position  errors  were  already  being  modeled  as  pseudo-range  errors,  then  this  shov'd  not  be 
overly  difficult. 
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figure  47:  heaoinc  error  autocorrelation  function 
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FIGURE  48:  PITCH  ERROR  AUTOCORRELATION  FUNCTION 
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FIGURE  49:  ROLL  ERROR  AUTOCORRELATION  FUNCTION 
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flQURg  SO:  HEAQINC/ROa  CROSS-CORRELATION  FUNCTION  fDETAILl 
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8.  CONCLUSIONS 


This  data  clearly  demonstrates  that  this  type  of  GPS  receiver  has  the  potential  to 
provide  very  useful  dynamic  attitude  information  at  a  reasonable  cost.  This  is  especially  true  in 
the  static  case  if  multipath  can  be  avoided  and  in  the  dynamic  case  if  inertial  aiding  is  available. 
The  data  collected  on  this  trial  demonstrates  that  this  receiver  can  provide  high  quality  attitude, 
at  a  one  Hz  rate,  with  80%  availability,  over  a  seven  day  period. 

The  statistics  for  the  discrete-dynamic  attitude  errors  (errors  at  the  3DF  one  second 
data  points)  are  given  in  Table  4,  repeated  here: 

TABLE  4.  DISCRETE-DYNAMIC  3DF  ATTITUDE  ERRORS 


SOURCE 

RMS 

ERROR 

(degrees) 

STANDARD 

DEVIATION 

(degrees) 

PERCE] 

TITLE  (degrees) 

raw 

68%  95% 

bias  removed 

68%  95%  100% 

HEADING 

0.47 

0.47 

0.08 

0.13 

0.04 

56.7 

PITCH 

1.14 

0.27 

1.13 

1.20 

0.05 

a 

ROLL 

1.28 

0.26 

1.27 

1.34 

0.05 

0.11 

6.5 

Perhaps  the  most  significant  statistics  here  are  the  95%  errors  with  bias  removed.  For 
the  heading  error  this  was  only  0.08  degrees  (1.4  milliradians)  and  for  the  pitch  and  roll  it  was 
0.11  degrees  (1.9  milliradians).  Although  the  bias  errors  were  much  more  significant,  they 
likely  could  have  been  removed  by  a  more  careful  installation/calibration  procedure.  The 
larger  RMS  and  standard  deviation  values  are  due  to  a  small  amount  (1%)  of  spurious  data. 

The  dominant  stochastic  error  model  parameters,  derived  from  the  valid  3DF  data  with 
bias  and  spurious  points  removed,  are  shown  in  Table  6,  repeated  here: 

TABLE  6.  DISCRETE  3DF  ERROR  MODEL  PARAMETERS 


MARKOV 

HEADING 

50 

0.1 

0.03 

0.022 

4,000 

PITCH 

5<|> 

1.0 

0.04 

0.035 

600 

ROLL 

8v|/ 

1.0 

0.04 

0.035 

600 

*  order  of  magnitude  only 
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There  were  also  various  other  smaller  components  identified  in  Section  6  above,  but 
the  most  important,  from  a  Kalman  filtering  point  of  view,  are  the  white  noise  and  Markov 
components.  These  important  components  are  all  about  the  same  magnitude,  between  0.022 
and  0.04  degrees  (0.4  and  0.7  milliradians).  The  correlation  time  for  the  pitch  and  roll  Markov 
errors  was  about  ten  minutes,  and  for  the  heading  was  about  one  hour.  Integration  with  an 
inertial  system  should  easily  remove  the  white  noise  component,  and  also  part  of  the  Markov 
components  (especially  for  pitch  and  roll,  which  have  short  correlation  times). 

There  were  some  significant  discontinuities  in  the  heading  error  corresponding  to  ship 
maneuvers,  suggesting  that  multipath  was  a  problem.  Also,  about  1%  of  the  data  was 
"spurious"  (errors  in  excess  of  0.2  degrees). 

These  observations  can  be  given  with  confidence,  since  they  are  such  that  the  errors  of 
the  inertially  based  reference  system  could  not  have  significantly  influenced  them.  This  is 
because  the  stochastic  behaviour  of  an  inertial  system  is  well  known  and  occupies  a  different 
part  of  the  autocorrelation  spectrum.  This  complementary  nature  of  the  error  characteristics  in 
fact  bodes  well  for  the  prospects  of  obtaining  a  synergistic  improvement  in  performance 
through  the  Kalman  filter  integration  of  these  two  types  of  system. 

It  should  also  be  mentioned  that  the  cost  of  this  type  of  receiver  is  expected  to  come 
down  in  the  future.  In  fact  other  models  of  this  type  of  attitude  measuring  GPS  receiver  are 
already  available,  at  substantially  lower  cost.  For  example  the  Trimble  "Tans  Vector"  is  being 
offered  for  less  than  $20K.  If  the  performance  of  these  lower  cost  receivers  is  comparable  to 
that  of  the  Ashtech  3DF,  as  determined  here,  then  they  will  make  very  attractive  sensors,  to 
complement  or  replace  inertially  based  systems. 

Of  course  significant  advances  are  also  being  made  in  fiber  optic  gyro  technology, 
which  could  improve  the  cost/performance  ratio  of  the  inertial  systems  as  well.  Since  the  two 
technologies  are  so  complementary,  this  leads  one  to  believe  that  more  combined  systems  will 
be  used  in  the  future.  There  are  already  various  inertial  systems  available  with  built-in  GPS 
positioning  receivers,  some  of  which  are  even  tightly  coupled  (optimally  integrated). 
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APPENDIX  A. 


PLATFORM  DYNAMICS 


Since  we  are  only  interested  in  the  angular  velocity  and  acceleration,  to  place  limits  on 
the  data  latency  and  interpolation  errors,  we  shall  assume  that  for  short  periods  of  time  the 
local  angular  behaviour  can  be  approximated  by  a  sinusoidal  variation  in  attitude  0,  (heading, 
pitch  or  roll)  presumably  induced  primarily  by  interactions  with  ocean  waves.  That  this  is  so 
can  be  clearly  seen  from  the  detailed  plots  of  Figures  12  to  14.  Therefore  we  shall  model  the 
attitude  angle  6  as: 

0  =  0O  +  i4sin(©/+<J>)  (Al) 

where  A  is  the  local  angular  amplitude  and  co  is  the  phase  rate  (in  radians/second).  Then  the 
angular  rate  dQ/dt  will  then  be 

—  =  co/4cos(co/  +  <j>)  (A2) 

dt 

and  the  acceleration  cPQ/dt2  will  be 

— 5-  =  -<n2.4sin(<D/  +  <j>)  (A3) 

dt 

Thus  we  can  place  limits  on  the  platform  dynamics  of: 

co  A 


®2A  (A5) 


We  can  see  from  the  data  (Figures  9  to  14,  for  example)  that  the  local  angular  amplitude.  A, 
during  this  trial  was  less  than  about  4  degrees  (and  largely  less  than  2  degrees): 

l  A  |  <  4  degrees  (A6) 

and  that  the  period  was  about  6  seconds.  Therefore  the  phase  rate  co  was  at  most  about 

co  =  2k  /6 

=  1 .05  radians/second  (A7) 

which  from  A4  and  A5  place  limits  on  the  angular  rate  and  acceleration: 
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and 


ae 

dt 


<  ( 1 .05  radians/sec)  (4  degrees) 
=  4.2  degrees/sec. 


<  (1.05  radians/sec)2  (4  degrees) 
=  4.4  degrees/sec.  ^ 


(A8) 


(A9) 
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APPENDIX  B. 


DAILY  ERROR  PLOTS 


This  appendix  simply  shows  the  discrete  attitude  errors,  with  biases  removed,  of  the 
valid  3DF  data  throughout  the  sea  trial.  This  is  the  same  data  as  shown  in  Figures  20  to  22, 
but  is  shown  here  in  24  hour  segments  so  that  the  temporal  behaviour  can  be  more  clearly 
seen. 

For  the  heading  errors  of  Figures  B1  to  B7,  a  second  plot  is  presented  in  which  the 
noise  has  been  suppressed,  using  a  simple  moving-window-averager,  to  more  clearly  illustrate 
the  Markov  component  and  the  discontinuities. 
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FIGURE  B1 :  DAY1  DISCRETE  HEADING  ERROR  VS  TIME 
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FIGURE  B2:  DAY2  DISCRETE  HEADING  ERROR  VS  TIME 
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FIGURE  B7:  DAY1  DISCRETE  PITCH  ERROR  VS  TIME 


FIGURE  B 1 1 :  DAY3  DISCRETE  PITCH  ERROR  VS  TIME 


FIGURE  B15:  DAY5  DISCRETE  PITCH  ERROR  VS  TIME 


FIGURE  B17:  DAY6-7  DISCRETE  PITCH  ERROR  VS  TIME 
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